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Abstract 

We study the asymptotic behaviour of the trace (the sum of the di- 
agonal parts) r n = T n (u>) of a plane partition ui of the positive integer n, 
assuming that u is chosen uniformly at random from the set of all such par- 
titions. We prove that (r n — con 2 ^ 3 ) /cin 1 ^ 3 log 1 / 2 n converges weakly, as 
n — > 00, to the standard normal distribution, where Co = £(2)/[2£(3)] 2,/3 , 
a = yVS/^S)] 1 / 3 and C(«) = V; . . 

DD 

1 Introduction 

Properties of various kinds of partitions are often studied using bivariate gen- 
erating functions of the following type: 

oo 

G(u,x;{a j } j > 1 ) = Y[(l-ux j )- a i =1+ ]T Q(m, n; Kb>i) u m x n . (1.1) 

j—l n,m>l 

Here u is finite, |ac| < 1 and {p,j}j>\ is a given sequence of non-negative numbers. 
A combinatorial interpretation of (|l.ip for integer-valued sequences {cij}j>i, is 
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obtained expanding the geometric progressions in the left-hand side. Setting 

oo a k 

A = U \J{k}^{X 1 ,X 2 ,...}, 

fc=U=i 

one can show that the coefficient Q(to, n; {a,j}j>i) equals the number of solu- 
tions in non-negative integers zj of the equations 

E x j z j = n, 
A36A 

< 

E zj = m. 
\ j 

Furthermore, a specific sequence of integers {dj}j>i may define a particular 
class of partitions, so that the nth coefficient Q(n; {aj}j>i) in the power series 
expansion of 

00 

G(l,x;{a j } j >i) = l[(l-xi)- a 3 = 1 + ^ Q(n; {a,-}^" (1.2) 

j=l n>l 

counts the number of partitions of n belonging to that class (see e.g. [U Chap. 
6]). A fundamental problem in the theory of partitions is to determine the 
asymptotic behaviour of Q(n; {cij}j>i) as n — > 00. In its most general setting 
this problem was studied by Meinardus [H] who introduced a set of assumptions 
on the sequence {ctj}j>i and obtained an expression for the leading term in the 
asymptotic expansion of this coefficient. Two important classes of partitions are 
covered by Meinardus' formula: (i) a,j = 1 and (ii) aj = j,j = 1,2,.... In the 
first case Q(n;{cij = l}j>i) equals the number of (linear) integer partitions of 
n (see (H Section 1.2]) and Meinardus' asymptotic formula implies the famous 
result of Hardy and Ramanujan [5] for the number of such partitions. In the 
second case Q(n] {aj = j}j>i) equals the number of plane partitions of n, whose 
asymptotic expression was previously obtained by Wright [27 . 

A basic restriction in Meinardus' scheme states that the Dirichlet's series 

00 

generated by the sequence {aj}j>i has to converge in the half-plane Sftes > a > 
0. It turns out that the asymptotic behaviour of the coefficients Q{m, n; {dj}j>i) 
in the bivariate power series expansion (jl.l[) strongly depends on the value of 
the parameter a. Haselgrove and Temperley [S] studied the case when a < 2, 
n — > 00 and m becomes large with n at a specific rate. They applied the classical 
method due to Hardy and Ramanujan [5] and obtained a result of the form 
of a local limit theorem for the ratio Q(m, n; {a,j}j>i)/Q(n; {o>j}j>i)- Their 
result established convergence to a non-Gaussian distribution. Haselgrove and 
Temperley [9l Section 3] also conjectured that the Gaussian law would appear 
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if the Dirichlet's series parameter a is not less than 2, however, a formal proof 
is still lacking. 

In this paper we consider the sequence cij = j,j = 1,2, . . ., and focused on 
plane partitions of positive integer n. A plane partition w of n is a representation 

n = H W M> 

i,3>l 

in which the array u> — (k>i,j)i,j>i of non- negative integer entries is such that 
u>i.j > Wi+ij and ojij > u>ij+i. We may also assume that to occupies the first 
quadrant of the coordinate system iOj. As an illustration the following is the 
plane partition of n = 8: 
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(1.4) 



O i 



For the sake of brevity the zeroes in (|1.4p are deleted, so that the abbreviation 

1 1 
3 2 1 

presents the plane partition of 8 in a shorter way. It seems that MacMahon [13] 
was the first who suggested the idea of a plane partition. Various properties of 
plane partitions and their applications to combinatorics, algebra and analysis 
of algorithms may be found in [U Chap. 11], [HI Chap. 12] and [22l Chap. 7]. 
Further, for the sake of brevity, we also let 

oo 

Q(n) - Q(n;{ aj =j}i>i), G{u,x) = [](1 - ux j yK (1.5) 

3=1 

In terms of notations (|1.5|) the generating function (|1.2|) becomes 

oo 

G(l,x) = U(l-^'r j = l + Y^Q{n)x n . (1.6) 

j=l n>l 

It is also easily seen that the sequence a 3 — j,j = 1,2,..., implies that 

oo 
3=1 

(see (|1.3p ). Therefore, in this case we have a = 2. 
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The main goal of our paper is to prove Haselgrove and Temperley's conjec- 
ture [9] Section 3] that the Gaussian law determines the asymptotic behaviour 
of Q(m, n) = Q(m, n; {a,j = j}j>i) if m becomes large with n at a specific rate. 

It turns out that this problem has an interesting probabilistic interpretation. 
If we introduce the uniform probability measure P — P n on the set of all plane 
partitions of n assuming that the probability 1/Q{n) is assigned to each plane 
partition w, then each conceivable numerical characteristic of uj becomes a ran- 
dom variable. An exact combinatorial expression for the numbers Q(n) does 
not exist. As it was previously mentioned their asymptotic was determined by 
Wright [57J and subsequently his result was confirmed by Meinardus' general 
theorem Q3]. It was shown that 



W - 2ll /3 6 3 3 i/2^i/2 ^ 25/36 exp{3[C(3)] 1 /3 (n/2)2 /3 + 2c }, (1 . 8) 

where 

r°° „,i™-„ 

dy. 



o 



3 2iry _ I 



(In the statement of Wright's main result [57J p. 179] the constant 3 1 / 2 in the 
denominator of (|1.8|) is missing, however, it is included in his final result at the 
end of the proof of his theorem on p. 189. Further results on the asymptotic 
expansion of Q(n) can be also found in [HE].) 

We will consider here the trace r„ of a partition cj defined as the sum of its 
diagonal parts: 

n 
3=1 

To study the asymptotic behaviour of r n as n — > oo we use the following gener- 
ating function identity (see [21] or [4] Chap. 11, Problem 5]) : 

oo n 



G{u, x) = 1 + Q(n) x 11 2^ P(r„ = m) u m = 

n— 1 rn— 1 

oo oo 

= 1 + Q( n ) Vn(«) x n = l[( 1 - «ar , ')" J '« (1-9) 



where 



ip. 



n {u)=Y J P ( T n=m)u m (1.10) 



is the probability generating function of the trace T n and G(u, x) is the generat- 
ing function defined by (jl.5p . It is now clear that if Haselgrove and Temperley's 
conjecture, that we stated above, is valid, then the trace t„ of a random plane 
partition has to be asymptotically normal as n — > oo. The main result of this 
paper is the following limit theorem. 
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Theorem 1 For any real and finite z, we have 



Inn P ( T " ~ ggff? < z ) = J= / e -. 2 /2 dyj 
™^°° VcinVsiog^^n / V27r 7 

— oo 

c = C(2)/ [2C(3)] 2/3 = .916597104, 
ci = v/T73/[2C(3)] 1/3 = .430977269. 

This result confirms Haselgrove and Temperley's conjecture when a = 2 and 
-D(s) = C(s — 1) ( see (jl.7[l ). We believe that our method of proof can be utilized 
to study the general case a > 2. 

Our study is also partially motivated by similar results obtained for linear 
partitions of n. Erdos and Lehner [5J were apparently the first who have studied 
random linear partitions using a probabilistic approach. As a matter of fact, 
they showed that the number of summands, after an appropriate normalization, 
converges weekly, as n — + oo, to a random variable having the extreme value 
distribution. Subsequent work by a number of authors provides considerable 
information about the structure of a "typical" linear partition of a large integer. 
(We refer the reader e.g. to [23], [2J], [25], [6], [7], [TS], [EE], [TSJ and the 
references therein.) To this effect, Theorem 1 continues the study of partitions 
initiated by the probabilistic approach. 

The proof of our main result is based on a method developed by Hayman 
[TO] . To get an estimate for the Cauchy integral stemming from (|1.9p we use 
the fact that the generating function G(l, x) defined by (|1.6p satisfies Hayman's 
admissibility conditions in a neighborhood of its main singularity x = 1 and 
outside it (see Lemmas 1 and 2, respectively). A relevant approach to problems 
related to other characteristics of random plane partitions may be found in [19] 
and [16]. 

We organize our paper as follows. Section 2 contains auxiliary facts on the 
admissibility of the generating function G(l,x) defined by (|1.6p and on the 
asymptotic behaviour of the coefficients Q{n) in its power series expansion. In 
Section 3 we present the proof of Theorem 1. 

2 Preliminary Asymptotics 

In order to get an asymptotical estimate for G(l, x) around its main singularity 
x = 1, we first introduce the analytic scheme of assumptions on the sequence 
{ a j}j>i of non-negative numbers, which is due to Meinardus [14] . The following 
three conditions must be satisfied: 

(Ml) The Dirichlet series D(s) (see (jl.3p ) converges in the half-plane 8e s > 
a > 0, and can be analytically continued into the half-plane 5Re s > — oto, 
ctQ € (0, 1). In 5ie s > — ao, D(s) is analytic except for a simple pole at s = a 
with residue A. 
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(M2) There exists an absolute constant a\ > such that D(s) = 0(\^sm s\ ai ) 
uniformly for 3?e s > — a>o as |Qto s\ — ► oo. 

oo 

(M3) Define g(v) = a,je~ JV , where i> = y + 2?:iw and y and w are real 

numbers. If |argu| > 7r/4 and |to| < 1/2, then 3?e g(v) — g(y) < —ct2y~ x for 
sufficiently small y, where x > is an arbitrary number, and «2 > is suitably 
chosen and may depend on \- 

We shall be concerned, in the first instance, with the behavior of G(l, x; {a,j}j>i) 
(see ll.2p as x becomes close to 1. Meinardus [13] (see also Andrews [U Lemma 
6.1] proved that under the assumptions (Ml) and (M2) 

G(l,e~ v ; {a^>i) = exp (^ir(a)C(a + l)v~ a - D(0) logv + D'(0) + 0(y a °)} 

(2-1) 

as y — > uniformly for | argu| < 7r/4 and |tu| < 1/2. (Here T(a) denotes Euler's 
gamma function and log (•) presents the main branch of the logarithmic function 
satisfying logv < for < v < 1.) 

Let us take now a sequence {r n } which, as n — > oo, satisfies 

r „ = 1 _M^ + M5^!_|fi + o ( „-n (2 . 2) 

n 1 '" 3 In^l 6 An 



For the sake of brevity we also set 



6C(3) 
(l-r)< 



b (r) = 7T^, (2.3) 



where < r < 1. It is easy to check that (I2.2[) and (|2.3p imply 

i, c»' = pfjF + 0( " ) (2 ' 4) 

as n — > oo. 

The next lemma suggests a tool that we shall subsequently use in Section 3 
to obtain the main term in our asymptotics. 



Lemma 1 If r n satisfies for large n, then 

G(l, r n e l9 )e- l9n = G(l, r^-^^l 2 [l + 0(1/ log 3 n)\ 
as n — > oo uniformly for \9\ < S n , where 

On = T 

logn 

and b(r n ) is determined by i2.3\) . 



(2.5) 



Proof. Our starting point here will be Meinardus' general asymptotic for- 
mula (|2.ip . We apply it for the sequence ctj = j, j = 1, 2, It is not difficult to 

show that A = 1 and -D(O) = —1/12. Classical result on the £ function implies 



G 



condition (M2) (see e.g. [551 Section 13.51 ]). Therefore, for v — y + 2iriw, we 
get 

G(l,e- v )=expU(3)v- 2 + ^ log v + D'(0) + O{y ao )^ (2.6) 
as y — > uniformly for \w\ < 1/2 and | &igv\ < 7r/4. Setting 



1 n c 7 



(2.7) 



we see that y — y n — ~^>gr n and io = —0/2tt. The asymptotic behaviour of 
— logr„ can be determined with aid of (|2.2|) as follows: 



logr„ 



[2 C(3)] 1/3 [2 C(3)] 2/3 C(3) 
[2C(3)] 2/3 2 4 / 3 [2 C(3)] 2/3 [C(3)] 1/3 2C(3) 



2n 2/3 



4/? 



3n 



+ 0(n- 4 / 3 ) 



[2 C(3)] 



1/3 



Combining 



n l/3 

, we observe that 



+ 0(n- 4 / 3 ), n 



G(l,r»e' e ) 
G(l,r„) ' 



(2.8) 



(2.9) 



y n - it 

Vn 



1/12 



exp {C(3)[(y„ - itf)" 2 - y,; 2 ] - *0n + O (y r ? )} 



A Taylor's formula expansion for |0| < S n yields 



(y n - i9)- 2 - y7 n 2 = 2i6y- 3 - 6-y^ + 0{\6\ 3 y^) 



2i9y- A - 30^" 4 + 0{8 A n y->) 



(2.10) 



Using (|2.8[) , we also get the following estimate for the factor outside the exponent 



~ I [2C(3)]V3 + 0(„-i) J 

= l + 0(5 n n 1/3 ). (2.11) 
Finally, we notice that (|2.8p implies the bound 

Vn° = 0(n- a ^ 3 ). (2.12) 

Hence, inserting j231), tpT8 |l - <[27l0jl - t[27T2]) into (|2T9]l . we obtain 

2i6»n 



G(l,r» e ' g ) _ l9n 
G(l,r„) 



l + O^n 1 / 3 ) ex P U(3) 



2C(3)[l + 0(n~ 1 )] 3 
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wm^mr^w + ° (fe5) J - l6n+0 ( n ~ W3 ) . 



l + 0(n- 2 / 9 /logn)l expK(3) 



iOn 

W) 



{l + 0{n- 1 )) 



~ 2^(3) [h{Tn) + ° (n)] ^ + 0(n_1) ) + °( <5 « n5/3 ) - i9n + ( n ~ a ° /3 
= 1 + 0{n~ 2/9 / logn)] exp {i(9n + 0(S n ) - 6> 2 6(r„)/2 + O (nSl) 
+0{8lb{r n )n- 1 ) + 0(^n 5 / 3 ) - i9n + O (V Q °/ 3 ) } 
= l + (9(n- 2 / 9 /logn)j exp{-6» 2 6(r„)/2 + 0(l/log 3 n)} 
= e -e%{r n )/2 + (1/ log 3 n)] . 

This completes the proof. 

We also need another lemma that will establish a uniform estimate for 
G(e zT ,x), \x\ — r n outside the range —8 n < argx < S n , if T is real and suitably 
bounded. 



Lemma 2 If r n and S n satisfy 
T = T(n) is such that 



and 



respectively, and the function 
(2.13) 



T = T(n) = e(n-^ 3 /y/logn) 
as n — > oo, then there exist two positive constants e and uq so that 

\G(e lT , r n e l( >)\ < G(l,r n ) exp [e - 2n 2 / 9 /[2C(3)] 4/3 log 2 n} 
uniformly for tt > \8\ > 5 n and n > n . 

Proof. By taking logarithm, for |x| < 1 and \u\ = 1, we get 

I X. ) oc 

k>gG(u, x) — log < — UX^)~^ > = — j log(l — MI 3 ) 



OO oo 



j=l 1=1 1=1 ] = 1 

Thus, substituting x — r n e %e and u = e lT , we obtain the estimate 



\G(e* T ,r n e w )\ = |exp {log G(e iT , r„e )} 



!oo 1 OO j 
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exp < 3 ' r n 



2=1 



Ey E-? 

1=2 j=l 



r U e il(jS+T) 



exp <^jV„ See 



exp 



oo oc j 

1=2 j=l J 

!OC OO 1 OO 

^ j r£ cos(j0 + T) + £ - ^ j r M co S [l(j6 + T)} 

2=1 i=2 2=1 



< exp J ^ 



j cos 



2=1 



where 



OO OO I 

(^+t)+Et e^' 

1=2 2=1 J 

{oo oo 1 OO 

5] j r> [cos(j0 + T) - 1] + ^ - £ j r« 
2 = 1 1=1 3 = 1 

= G(l,r„)exp T)}, 



H n (9,T)=Y,jri [cos(j0 + T)-l] 

3 = 1 



5ie 



i(9+T) 



(1 - r n e iS f 



(l-r„) 2 



_ r„ [cos(0 + T) + r,2 cos(6> — T) — 2r„ cos(T)] 
~ (1 - 2r„ cos(6>) + r2 )2 

It is not difficult to show that the function 



(1 - r n f 



(2.14) 



(2.15) 



f(9) = cos(6> + T) + r£ cos(6> - T) 

in the numerator of H n {6, T) attains its maximum value in the range 8 n <| 9 |< 
7r at = ±<5„. To prove this one needs to determine the behaviour of f(9) in 
the ranges T < 6 < tt and S n < 6 < T. The first case is an easy exercise whose 
study avoids the asymptotic form (|2.13|) of T. In the second case one can use 
the following representation of the derivative 

f'{9) - - sin(0 + T)-r 2 n sin(<? — T) = -(1 + r 2 n )9 - (1 - r 2 n )T + 0(T 3 ). 

dH2J) and (|2~T3|) show that 

(1 - r 2 n )T = 9(7i- 2/3 log- 1/2 n), 0(T 3 ) = O^ 1 log~ 3/2 n). 



9 



Therefore, for sufficiently large n, the value of f'(9) becomes close to — (l+r%)9. 
This implies that f(9) decreases for 8 n < 9 < T and increases for —T < 9 < —S n . 
It follows that one may replace 9 by — <5„ in (|2.15[) . Thus we can write 



,T) < 



r n [cos(-6 n + T) + r 2 cos(-£» — T) — 2r n cos(T)] 



(l-2r n cos(<y + r2)2 



(1- 



(2-16) 

To estimate the cosine functions in the right-hand side of this inequality we 
express 5 n and T by (|2.5[) and (|2.13p , respectively We get the following expan- 
sions: 

cos§n = i_^ + O ((5Ji) = l-^- 3 _ + O(n- 20 / 9 log- 4 n), 
2 2 log n 



cosT± 2sm(T=F -y) sm "2"' 



cos(-5 n ± T) 

cosT = 1 - e(n~ 2/3 /logn). 
Moreover, (12. 2|) implies that 

(l-r„)~ 2 = 



2C(3) 



2/3 



+ 0(n 1 / 3 ). 



Substituting these estimates in the right-hand side of (|2.16p , after some manip- 
ulations, we obtain 



H n (9,T) < 



r n [l-e(n-2/3/l gn)] 



+ (l-r„) 2 log 2 n ~i~ u \ (1-r, 



) 2 log 4 i 



r n [l-6(n- 2 / 3 /log«)] 



(1 - r„) 2 [1 + 2n- 4 /9/[2C(3)]2/3 log 2 n + 0(n-W / log 2 »)] (1 - r„) 2 



r ra [l-e(n- 2 / 3 /logn)] f 2n~ 4 / 9 
(1 - r„) 2 ' - 



2C(3) 



2/3 



[2C(3)]2/3 1 g 2 n 

2n- 4 / 9 



O 



log 2 n 



2n 2 / 9 



[2C(3)] 2 / 3 log 2 n 
(9(1/ log n). 



(l-r„) 2 
(9(1/ log n) 



[2C(3)]4/3log 2 n 

Inserting this into (|2.14p . we obtain the required bound. 

Further, we shall essentially use the asymptotic form of the numbers Q(n). 
It is given by Wright's formula (|1.8[) . however, we need this result in a slightly 
different form. It is not difficult to verify that Lemmas 1 and 2, (|2.2p and (|2.3p 
imply that G(l,x) belongs to the class of Hayman's admissible functions [10] 
and thus, Hayman's general asymptotic formula for the coefficients in the power 
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series representations of admissible functions is valid for the number of plane 
partitions of n, Q(n), as well. The next lemma encompasses these results. We 
only sketch its proof and insert a remark explaining the role of the asymptotic 
expansion (|2.2p there. 



Lemma 3 For G(l,x), the generating function of the numbers Q{n) of plane 
partitions of n, defined by U.6\) , we have 

Q(n) ~ G(l, r„)r-"/[2^6(r„)] 1/2 (2.17) 

as n — > oo 7 where r n satisfies the equation 

rG'{l,r)/G{l,r) =n (2.18) 

for sufficiently large n and b(r n ) is defined by 



Sketch of the proof. It is clear that \x\ = 1 is a natural boundary for G(l, x). 
Lemma 1 shows the behaviour of G(l,x) around its main singularity x = 1 
(condition I of Hayman's Definition [TU]); Lemma 2 establishes the negligibility 
of the growth of G(l,x) as x — ► xo, \xo\ — 1 and x n ^ 1 (condition II of [TU]). 
It is then easily seen that 

rg(l,r)_ 2 ^ r* h g r? 



G(l,r) ^(l-ri)3 ^ (1 - r^) 2 
+ ((l-r)" 3 ) 



2C(3) , _„., ^_ 3 , 



(1-r) 3 

as r -> r. This enables one to conclude that r„, determined for sufficiently 
large n by (|2.18p . can be substituted by the asymptotic expansion (|2.2p . Thus 
one can obtain (|2.17|) after a direct application of Hayman's theorem [10] . 

Finally, we notice that (|2.2p - (|2.4j) and (|2.6p imply the coincidence of the 
right-hand sides of P]) and ([2~T7]) . 

3 Proof of the Main Result 

First, we let in PI) 



G(u, i) = e^*), (3.1) 

that is, we set 

oo 

F(u,x) = -Y / j l og(l-ux j ). (3.2) 

i=i 

We now apply Cauchy's coefficient formula to (|1.9p on the circle x = r n e %e , 
— 7r < < 7r, with r„ determined by (|2.2p . Thus, in terms of the notations 
(fTTTUj) . dHUJ) and (021), for |u| < 1, we obtain 

Q(n)v?„(u) = - s — exp {F(u, r„e ie ) - i6»n} d6>. 



li 



We break up the range of integration as follows: 

Q{n)ip n (u) = Ji(n, u) + J 2 (n, u), 

where 



Ji(n, u) 



J 2 (n, u) 



2tt 



-5„ 



exp {F(u,r n e id ) - iOn) d0, 



2tt 



/ exp {F(u, r n e i0 ) - iOn) dO. 

JS n <\6\<TT 



(3.3) 

(3.4) 
(3.5) 



3.1 An asymptotic estimate for Ji(n, u). 
Using Taylor's formula expansion, we can write 

F(u, r n e w ) = F(u, r„) + r n (e* e - 1) ^-F(u, x) 



X TV 



f(e*-lf —F(u,x) 



9 3 



(3.6) 



To evaluate the partial derivatives of F(u, x), we shall use the following expres- 
sions: 

d 



dx FM 



. 1 1 - ur J n 



x = r 



n 



dx 2 



F{u,x) 



i 2 (i 1W- 2 „-3 20-1) 
: h M > — 

._ , i - uri £i a 



J=2 



j 3 (j -l)(j-2)rJ~ 3 
1 - «ri 

^ ,-4 3(i-l) 



9a; 



oo 



J=3 



f^ 2 (i - f-r (1 



(3.7) 
(3.8) 



3 t~L C 1 ~ Mr «) 3 

We proceed further to establish the convergence of Ji(n, u) in terms of 
Fourier transforms by setting u — e lT , — oo < T < oo, in (|3.4[) . (|3.6p - (|3.9|) . 
In what follows later, an application of Levy's continuity theorem for charac- 
teristic functions [12j Section 3.6] would specify the value of T as a function of 
the main parameter n. Furthermore, we need notations for the following four 
functions: 



(3.9) 



( e y-* T - i)> 

(k,m) = (0,1), (1,2), (1,3), (2, 3), < z < oo. 



(3.10) 



They are closely related to the Debye functions (see e.g. [T| Section 27.1]). 
Moreover, formula 27.1.3 of .1! shows that 



VVm-i(0,0) =m!C(m + l). 



(3.11) 
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We now proceed to the asymptotic estimates of the summands in 
Interpreting again the sum in (|3.7|) by a Riemann's one with the same step size 
y n = — log r n as in the proof of Lemma 2 of [TS] and using (|3.1ip and (|2.8[) , we 
find that 



r n ^F(u,x) 



x = r n ,u = e" 
[^,i(0,T) + O(n- 2 / 3 )]y- 3 + 0(^/3) 
: [^2,1 (0,0) + R(n,T)]y- 3 + O^ 1 / 3 ) 



= [2C(3) + R(n,T)} 



2C(3) 



■0(1) 



(3.12) 



n+^^ + 0(\R(n,T)\) + 0(n^) 



where 



R(n,T)=T^L 1 (0,T 1 ) = iT 



00 y 2 e y- iTl dy 



-iTi 



I) 2 



is the remainder term in the Taylor's formula expansion of ^2,i(0,T) and < 
T\ < T. Since the last integral is finite, from (|2.13|) we get 



R(n,T) = 0(n- 1 / 3 / v 4oi^). 



(3.13) 



Furthermore, note that (|3.6[) requires a multiplication by e l6 — 1. Hence by 
(I2~2l . (l2~5l) (I3~T21 (E33D 



r„(e l9 -l) ^f(^) 



a; = r„ , m = e 



(3.14) 



^ + O(|0| 2 )] 



nR(n,T) m , 

" + ^3r +0(n 



= + i 9 n f*P + 0(S n n 1 / 3 ) + 0(n<5 2 ) 

= zfa + »g n ^ 3 j r) +0(n- 1 /Vlog 2 n). 

To deal with the third term of the expansion in (|3 . 6|) , one has to follow the same 
line of reasoning. Thus, in a similar way (|3 . 8[) becomes 



dx 2 



F(u,x) 



yn 4 [MVn,T) + i> 5 ,i{Vn, T)] + 0(n). (3.15) 



■x—r n ,u—e z 



On the other side, using Taylor's formula expansion and (|3.1ip . we obtain 



VMl/n,T)=lM0,T) + O l/ n 



oz 



(3.16) 
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^3, 2 (0, T)+0(l/n) = ^3, a (0,0)+O(|T|)+O(l/n) = eC^+Ofa-^/Vlog^), 



^3,l(2/n,7 1 )=^3,l(0^) + O^„ 

= ^3,1(0, T) + 0(n- 2 / 3 ) = ^3,1 (0,0) + 0(|T|) + 0(n- 2 / 3 ) 



(3.17) 



= 6[C(3)-C(4)]+0(n- 1/3 /v / l^)- 
So Q3) and (l3"T5D - (IBTTD imply 



(3.18) 



= l{l + 0(n-V 3 )][-9 2 + 0(\9\ 3 )} 



x—r n ,u—e zT 

4/3 



/ n 



0(n^ 3 ) 



V2C(3) 

x{6C(4) + 6[C(3) - C(4)] + 0(n- x / 3 /v^i^)} 

1 r / \ 4 / 3 

= i[l + 0(n- 1 /3 )]M 2 + 0((5 3 )] +0(n i/ 3) 



-02 / „ 
T~ ^2C(3) 



4/3 



x[6C(3) + (9(n- 1/a /Vlog^)] 
[6C(3)][1 + Oin- 1 ' 3 )] + 0(n 1/3 <5 2 ) + 0{n 4 ' 3 8 3 n ) 



Q2 / \ 4 /3 



2 V 2 C(3) 
? 2 / n \ 4/3 



[6<(3)] + 0(n 4 / 3 S 3 ) + 0(nS 2 n ) 



Uc(3)J 



Finally, following similar but simpler analysis as above, with the aid of (|3.7p we 
can show the negligibility of the error term in (|3.6|l . We have 



/ d 3 
0(\9\ 3 g^F(\u\, X ) 



0{5 A n y. 



3„,-5\ 



0((n- 15/9 /log 3 n)n 5/3 ) = (9(1/ log 3 n) 



(3.19) 



We are now ready to substitute (|3J]1 . (j3~T4)l . ([3~T8]) and (f3~19|) into the 
integral of Ji(n, e lT ) (see the p.4[l ). We can write 



4/3 



2tt 



«p|<fln + itf- 2C(3) 



2((3) 



[6C(3)] +0(n- 1 / 9 /log 2 n) +0(1/ log 3 n) - »0n 
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F(e 



27r6 1 /2( r „) 



1 + 



1 



log 3 n 



5 " fel/2(r " ) f itn 1 / 3 R{n,T) f 



exp 



[2C(3)]5/6 v / 3 



'-5„b 1 /2( r „) 

(3.20) 

Note that we substituted 9 — t/6 1 / 2 (r„), — oo < t < oo, with 6 1 ^ 2 (r„) defined 
by (|2T3|) and (|2~4|) . In addition, (|2~4| and p3| justify the computation 



(3.21) 



So, the bounds of the last integral in (|3.20|) tend to ±oo. The estimate (|3 . 1 3|) for 
R(n, T) and the Lebesgue's dominated convergence theorem allow the passage 
to the limit under the integral. Moreover, (|3.13p shows that the limit of the 
integrand equals e - * 2 / 2 with an error term of order 0(l/\/logn). The addi- 
tional error term that we get replacing ±<5„6 1 / 2 (r„) by ±oo can be estimated 
using the asymptotic expansion of the function (27r) -1 / 2 J" z °° e _t dt as z — ► oo 
(see [TJ Chap. 7]). Therefore, (|3.2ip implies that this error term is at most 
0(n _1 / 9 (log n) exp{— d 2 n 2 / 9 log 2 n}) = 0(l/y/Togn). This shows in turn that 



Ji(n,e iT ) 



-n F(e' J ,r„) 
1 n c 

27T6 1 /2( r „) 

r -n e F(e iT ,r n ) 



l + 0(l/v1ogn) 



~ dt 



l + 0(l/Vlogn) 



(3.22) 



V2^«) 

This completes the required estimate for J 1 {n,e lT ). 
3.2 An asymptotic estimate for J 2 (n,u). 

To estimate this integral asymptotically we only need to apply directly Lemmas 
2 and 3. So for n > no and bounded T, we get from p. 41) and (|2.4j) : 



\J 2 (n,u)\ < 



2n 



G(e lT ,r n e l0 )\d9 



>\e\>8„ 



< r- n G(l, r n )(l + e) exp{-2n 2 / 9 /[2C(3)] 4/3 log 2 n} 
r- n G(l,r n ) i 



(3.23) 



■{1 + e) exp{-2n 2 / 9 /[2C(3)] 4 / 3 log 2 n + log ^2Trb(r n )} 



~ 0»(1 + e) exp{-2n 2 / 9 /[2C(3)] 4 / 3 log 2 n + - log ? i + 0(1)} = o(Q(n)). 

3.3 Final estimate for Ji(n, e lT ): the choice of T = T(ri) 
We start our estimation with a Taylor's formula expansion: 



d 



F(e 11 , r n ) = F(l, r n ) + (e 11 - 1) —F(u, r. 
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+(e lT ~l) 2 |^(«,rn) 



u=l 



O \T\ C 



Q3 



u=l 



F(l,r n )+iT —F{u,r n ) 
ou 



T 2 ^F(u,r n) 



(3.24) 



+0 \T\ 



a 3 



d 2 



+ \Tf 



The partial derivatives of F(u, r n ) can be again evaluated by Riemann's integral 
in a way similar to that presented in the proof of Lemma 2 of [15] . We have by 
D, (j3TTU|) and (|3~TT|) that 



d_ 

du 



u=l i= i 



Tig 



2C(3) 



2/3 



[l + CKn- 1 )] Vi,o(0,0) l + 0(«- 1/3 ) +0( ? i 1/3 ) 

2/3 

C(2) + ©(n 1 / 3 ). 



2C(3) 



Similarly 



a 2 



E 



2j 



= Vn 2l {Vn) +o{yll(y n )), 



(3.25) 



(3.26) 



where 



-2r 



dv = 



!Jn 



dv 



l yn (1-e-) 2 "' e»--l 
Since the first two summands above have finite limits and 

f°° dv 

Jy n e»-l =1 + q( 



'(1) 



as y n — > 0, we observe that 

% n ) = 0(1) - logy n [l + o(l)] = [1 + o(l)](- logy„). 



Substituting this into ()3.26j) and invoking the asymptotic (|2.8|) for ?/„, we find 
that 



du 2 



F(u,r n ) 





n 


u=l 





-i 2/3 



logn 



o (n 2 / 3 log 



(3.27) 



In the same way one can deduce the following estimate for the third partial 
derivative: 

d 3 



du 3 



F{u, r n 



O(n). 



(3.28) 
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Putting (pT25j) . (j3~27) and (pHS)) into we obtain 



F(e lT ,r n ) = F(l,r„)+ir 



2C(3) 



2/3 



C(2)+0(|T|)- 



T 2 



2C(3) 



2/3 



o (T 2 n 2 / 3 logn) + 0(n|T| 3 ) + 0(n 2 / 3 T 2 ). 



logn 
3 

(3.29) 



It is now clear how to choose the variable T = Tin) which satisfies (12 . 13[) . 
Note that this choice determines the scaling factor in our limit theorem. Setting 



T = w 



2C(3) 



1/3 



logn 



-00 < W < CO, 



(3.30) 



(in agreement with (|2.13p ) we can rewrite (|3.29p in the following form 



F exp < iw 



F(l,r n ) + iw 



2C(3) 

-,1/3 



2C(3) 



C(2) 



logn 



log n w 



(3.31) 



o(l) 



The convergence here is uniform with respect to w belonging to any finite inter- 
val. Going back to the estimate Q3.22p for Ji(n,e iT ) and replacing there (|3.30p 
and (|3.3ip . we find that 



Ji n, exp < iw 



-l 1/3 



2C(3) 



l + 0( 



=) 



Vlogn 
1 Q(n) exp 





n 


exp | iw 


L2C(3)J 



logn 
3 

1/3 



-n p F(l,r n ) 



y/27rbU 



(3.32) 



C(2) 





n 




[2C(3)J 



-1 1/3 



C(2) 



log n u; 2 
~~3 2 



where the reduction in the last asymptotic equivalence is due to a direct ap- 
plication of Lemma 3. To handle with J\ in a more appropriate way we may 
rewrite (|3 . 32[) as follows: 



Ji yn, exp <^iw 

1/3 



2C(3) 



1/3 



logn 



x exp < —iw 



2C(3) 



C(2) 



logn 



Q(n)e~ 



2 /2 



(3.33) 
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It remains now to deal with the characteristic function of (t„— con 2 / 3 )/ cin 1 / 3 log 1 / 2 n. 
The expressions for cq and c\ in the Theorem and (|1.9p imply that it is equal to 



exp 



2C(3) 



1/3 



C(2) 



Xip n exp < !U) 



2C(3) 



1/3 



logn 



logn 



Setting u 



exp 



c n 



1/3 



1 1/2 

ci log ' n 



</5„(exp{iw/cin 1/3 log 1/2 n}) = $„(«;) 



in (|3.3p with T given by (|3.30|) and substituting estimates (|3.32l) 



and (|3.23|) in it, we obtain 

Q(n)#„(«0 = Q(n)e^ 2 / 2 + o{Q{n)). 

The required week convergence follows from Levy's continuity theorem for char- 
acteristic functions p~2j Section 3.6]. 
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